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Using exact optimization methods, we find all of the ground states of (±/i) random-field Ising 
magnets (RFIM) and of dilute antiferromagnets in a field (DAFF). The degenerate ground states 
are usually composed of isolated clusters (two-level systems) embedded in a frozen background. 
We calculate the paramagnetic response (sublattice response) and the ground state entropy for the 
RFIM (DAFF) due to these clusters. In both two and three dimensions there is a broad regime 
in which these quantities are strictly positive, even at irrational values of h/J (J is the exchange 
constant). 



Disordered magnets provide a paradigm for disor- 
dered systems in general, and they continue to be inten- 
sively analyzed by a variety of methods. However, due to 
metastability, conventional (e.g. Monte Carlo) analysis 
Q of the equilibrium domain structure of random mag- 
nets is often unreliable, especially at low temperatures. 
Since the ground state behavior is an important indicator 
of the low temperature behavior of most random magnets 
||, exact methods for ground state analysis are desirable. 
Fortunately, the true ground states of random magnets 
can often be found using optimization methods. 

The relation between optimization and random mag- 
nets was pointed out some time ago ||,[|. However ex- 
tensive use of these methods is more recent, partially due 
to the availability of more efficient algorithms. An exact 
optimization procedure to find the random field ground 
state was implemented by Ogielski ||. More extensive 
analyses on larger system sizes have been published re- 
cently These methods have also been extended to 
the analysis of the ground state degeneracy of random 
magnets [|]||]. Here we present, using a new algorithm, 
a more precise analysis of the ground state degeneracy 
and its consequences in RFIM and DAFF in dimensions 
d = 2,3 (square and cubic lattices). We concentrate on 
the following three aspects of of these degenerate ran- 
dom magnets: (1) The degenerate domain structure of 
RFIM ground states (e.g. Fig. 1). (2) The order param- 
eter which couples to the ground state degeneracy. (3) 
The ground state entropy, and in particular the physical 
origin of its continuous and discontinuous parts. 

We consider the random field Ising model (RFIM) with 
a binary random field, 

Hrfim = - J ^ OiV] — ^ Kai (1) 

where hi — ±/i, with h and J positive and the plus and 
minus random fields occur with equal probability. We 
also analyze the dilute antiferromagnet in a field (DAFF), 

Hdaff = / ] JxiXjajaj - hy XjQ-j (2) 
(y) i 

where pi = pS(xi — 1) + (1 — p)5(xi) is the probability that 
a site is present. In both cases, we analyze the ground 



state properties as a function of the ratio H = h/J on 
square and cubic lattices. In the DAFF case there is the 
additional parameter p which we fix at p — 0.9. 

The ground state critical behavior of random field mag- 
nets is still not completely understood [p|PJlO[|. At small 
H, large ferromagnetic domains are favored, while at 
large H > H*, the spins freeze along the directions of 
the local field hi. In two dimensions there is no sponta- 
neous magnetization for any finite H(H 2d = 0), though 
there is a rapidly growing ferromagnetic domain size 
I ~ exp(l/H 2 ) which can masquerade as a phase tran- 
sition at H <~ 1. In three dimensions, there is a spon- 
taneously magnetized state at small H < H c . Although 
the field theory analysis and early simulations suggested a 
continuous behavior in magnetization m(H) as H — > H~ , 
precise numerical work using exact optimization methods 
suggests a large jump in m at H c (for the (±/i) random 
field case H* d = 2.21 ± 0.01 and Am - 0.8) j§. The 
DAFF was introduced as a possible experimental real- 
ization of the RFIM |ll| , and an extensive literature has 
developed from this observation The DAFF is an 

antiferromagnet at small H, and at large H > H* all 
spins are polarized in the field direction. The DAFF sub- 
lattice (staggered) magnetization is believed to be qualita- 
tively similar to that of the magnetization of the RFIM. 
Our calculations are for the ground state degeneracy in 
the non-trivial regime < H < , where H* = 2d (d is 
dimension) is the field amplitude beyond which all spins 
follow the local field direction and the ground state is 
non-degenerate (for both RFIM and DAFF). 

The ground state degeneracy of the RFIM has been 
intensively studied in one dimension jl3) . There has also 
been an analysis on Cayley trees |L4| , and an interesting 
analysis on the square lattice The latter paper did 
not have the advantage of exact optimization methods 
and missed some of the key features of the ground state 
degeneracy. Although the infinite range model misses 
entirely the degeneracy we find here Jl4| , the one dimen- 
sional and Cayley tree models have several qualitative 
similarities with our results. More recently Hartmann j^] 
has presented a low precision calculation of the ground 
state degeneracy of random magnets, though the physics 
we elucidate here was not discussed by him. 
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In order to find the ground state of RFIM and DAFF, 
we use the mapping of these problems to a flow problem 
in combinatorial optimization (so called min-cut/max 
flow) Jl(|. This algorithm also gives the exact minimal 
energy surface in random networks 0,^8| ■ This has been 
known for some time however improved algorithms 
(push-relabel with global updates [^9|) now allow opti- 
mization of 100 3 lattices in a few minutes on a high end 
workstation. Our method relies on the concept of resid- 
ual graph introduced by the network flow algorithms [fl9[ . 
The full residual graph of the equivalent network flow 
problem holds the whole information about the ground 
state structure. A naive search of the ground states - 
which is equivalent to finding the domains that can be 
flipped without altering the energy (or the max-flow in 
the network flow terminology) - is exponential. Instead 
we generate a supergraph which shows how the domains 
are related to each other, i.e. which domains can be 
flipped independently. It turns out that many of the 
domains are independent and the exponential search is 
reduced to the few remaining dependent domains, and 
we search over these remaining domains. We show that 
the problem is equivalent to finding all the directed cuts 
in a directed graph with single arcs and no cycles. A 
typical supergraph is shown in Fig. [j]. It is easy to see 
how the structure of the supergraph makes such a search 
effective, as it reduces it to searches over much smaller 
independent graphs. Details of the method will be pre- 
sented elsewhere p0[ . 

A typical ground-state domain structure of the two- 
dimensional ±/i RFIM is presented in Figs. ||. Here 
green domains are fixed in the direction of the positive 
fields (dots), while white domains are fixed in the oppo- 
site direction. Domains of any other color can be flipped 
without changing the ground state energy (note that not 
all can be flipped independently, but the dependent do- 
mains are organized in clusters that can also be flipped as 
a whole). These domains produce a finite ground state 
entropy. Surprisingly, domains that can be flipped in 
the ground state exist even for irrational H — h/J. 
Note that these domains occur at the interfaces between 
the up-spin and down-spin domains of the RFIM ground 
state. The degenerate clusters at irrational H have zero 
field energy and the same exchange energy in both the 
up and down states of the cluster. For the RFIM on the 
square lattice, the lowest order degenerate clusters of this 
sort are indicated in the inset to Fig. |^. The number of 
these clusters (or two- level-systems (TLS)), utls, can be 
estimated using simple arguments: 



T-cut 



n T LS oc ptls — 



(3) 



where L is the linear size of the system, I oc exp[(l/H) 2 ] 
is the typical size of the ordered domains || and L 2 jl is 
the total length of interface between up and down spin 
domains in the system. Ptls is the probability of occur- 
rence of a TLS at a given interface site, ptls = Pn/4, 




S-cut 



FIG. 1. Typical supergraph as obtained using the algo- 
rithm. S is the set of spins frozen up and T the set of spins 
frozen down. Ci, C2 and C3 are independent clusters; they are 
made of subclusters that are not independent of each other. 
The S — cut is the ground state with all the independent clus- 
ters down, and the T — cut is the one with all of them up. We 
also show a directed cut (ground state) in which part of each 
independent cluster is up/down. 

where 1 /4 is the probability of occurrence of the up- 
down pair of fields and p n is the probability to have this 
pair surrounded in the ground state by frozen spins with 
the appropriate configurations. The entropy density is 
then s cx p n exp[- (1 / H 2 )] for H < 4 and for H > 4. p n 
is discontinuous at H = 2, because the dominant TLS's 
for H < 2 arc different than those for H > 2 (for example 
there are twice as many spin configurations that lead to 
a TLS below than above H — 2). If we take the observed 
jump ~ 3.3, then the above argument leads to the curve 
given in the inset to Fig. [3, which is very close in form 
to the continuous part of the entropy presented in Fig. 

I 

The series of sharp peaks occurring at rational values 
of H are due to additional degeneracy occurring when 
clusters have same value for the field energy plus the 
exchange energy in either the up or down states. These 
peaks can only occur at rational values of H and the 
cluster geometries which contribute at each rational are 
different. Naturally high order rationals correspond to 
complex clusters and have greatly reduced degeneracy. 
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FIG. 2. Typical ground states of the d — 2 random field 
ising model (Eq. 1) with H = h/J = 3/2 (top), H = 2 
(middle) H — 5/2 (bottom) (system size 50 x 50). Green 
indicates an up spin, white a down spin, and the other colors 
indicate spin clusters that can be flipped up or down (but not 
all independently of each other) without changing the ground 
state energy. Black dots indicate the sites where the random 
field favors the up spin orientation. 

It has been suggested before jl5| that for the 2d RFIM, 
the highest degeneracies occur at rationals H n = 2 + 2/ 'n 
for 2 < H < 4, with n = 1, 2.... Using our algorithm we 
checked this idea by considering n = 1 , . . 1 1 and all ra- 
tionals with denominators 2, 3, 4, 5, 6. We find that those 
with H n = 2 + 2/n are indeed dominant |2l|], see Fig. [^, 
and there is a similar sequence at fields 4 + 2/n in 3d (the 
zoology of the clusters leading to the dominant peaks is 
straightforward though tedious to enumerate). In the 
regime < H < 2 the 2-d RFIM ground-state entropy 
has spikes at a large number of rationals (see Fig. ||). 
These features are quite similar to those found in one di- 
mension (see Fig. 4 of [^3|) and on Cayley trees (see Fig. 
4 of jl4|] ). We have also done a preliminary analysis of 
the 3-d RFIM and the 2-d and 3-d DAFF ground states. 
In general we find that the RFIM and DAFF magnets 
in 2-d and 3-d are massively degenerate in the regime 
H c < H < and that their ground state entropy is 
finite even at irrational H . 
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FIG. 3. The ground state entropy of the 2-d random field 
Ising model (RFIM) as a function of h/J — H. The inset 
shows the smallest two-level-systems (TLS's) at irrational H 
and an estimate of the entropy (from Eq. (3)) produced by 
them. "+" indicates an up spin "-" indicates a down spin and 
a dot indicates where the local random field favors the "+" 
spin direction. The system sizes used were from 10 x 10 to 
130 x 130 and the entropy was found as the slope of the line 
(InD).vs.N, where D is the degeneracy, N is the system size 
(total number of spins) and the average is over the disorder 
(1000 samples were used). 

In the regime H c < H < H*, we can consider the 
ground state to be composed of a frozen background in 
which is embedded a set of largely non-interacting free 
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superspins (corresponding to each independent cluster). 
The ground state of these magnets thus can be considered 
to contain a large number of magnetic two-level-systems 
p2fl . There is a paramagnetic response at low tempera- 
tures for both the RFIM and DAFF in this regime. The 
natural ground state order parameter for the paramag- 
netic response in the regime H c < H < H* is the magne- 
tization mpferro for the RFIM and the staggered magne- 
tization mpstaggered for the DAFF. It is straightforward 
to calculate these order parameters using the exact op- 
timization algorithm, either by applying an appropriate 
infinitesimal field or by polarizing all of the degenerate 
domains in a given orientation (we do the latter). The 
results for the RFIM are presented in Fig. ^ for both 
square and cubic (inset) lattices. 




FIG. 4. The order parameter for the ground state paramag- 
netism (mpf erro ) of the RFIM on square and cubic(inset) lat- 
tices. System size was 200 x 200 for the square and 40 x 40 x 40 
to 60 x 60 x 60 for the cubic lattice (1000 samples were used). 

It is seen that the basic features of the ground state de- 
generacy (Fig. |^) are reflected in the ground state param- 
agnetic magnetization. In 3-d, the entropy remains zero 
at low H (at least for H irrational), reflecting the ferro- 
magnetic state for H < H c ~ 2.21 §. For H > H* = 6 
(in 3-d) , the spins are aligned with the random field, and 
the ground state is non-degenerate. Note that in addi- 
tion to the paramagnetic magnetization, there is a spon- 
taneous magnetization msf erro for H < H c . In experi- 
ments in which the ferromagnetic field is swept to pro- 
duce a magnetization loop, the measured zero field mag- 
netization is the sum i.e. m0f erro = msf erro + mpf erro . 
Thus there is a finite equilibrium magnetization jump 
at zero temperature even for H > H c . Of course 
mpferro = for T > 0, but the effects of the ground state 
degeneracy should be reflected in magnetization anoma- 
lies and a Curie Law in the susceptibility at low temper- 
atures. 

Calculations of the paramagnetic order parameter for 
the DAFF, mp stagge.re.d-, for square and cubic lattices is 
presented in Fig. pL 




FIG. 5. The order parameter for the ground state paramag- 
netism (mpstaggered) of the DAFF on square and cubic (inset) 
lattices. The same system sizes and number of samples as for 
the RFIM were used. 

Qualitatively the situation is similar to that in the 
RFIM. There is a strong sublattice paramagnetic re- 
sponse for all H c < H < if*, with spikes at certain 
rational values. These figures are for a DAFF with dilu- 
tion p = 0.9, but only the details change as p is varied. 
In the regime H < H c there is a spontaneous staggered 
magnetization, and low temperature measurements (such 
as neutron scattering and NMR) should be influenced 
by both the "staggered paramagnetic" response and the 
spontaneous staggered magnetization. We also note that 
the existence of the additional order parameter mpf erro 
in the case of the ±h RFIM suggests that the ±h RFIM 
and the Gaussian RFIM may not be in the same univer- 
sality class. 

We have described two developments in the analysis of 
random magnets. 

(i) Using optimization methods it is possible to efficiently 
calculate the ground state structure of RFIM and DAFF 
magnets (see Fig. 2). (ii) The (±h) RFIM and the 
DAFF magnets have a spontaneously ordered state for 
H < H c , a massively degenerate ground state in the 
regime H c < H < and a non-degenerate ground state 
for H > if* . In the degenerate regime there is a strictly 
positive paramagnetic response and ground state entropy 
even at irrational H, with additional degeneracies at ra- 
tional H (see Figs. 2-4). 
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